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We discuss a general framework for the realization of a family of Abelian lattice gauge theories, i.e., link 
models or gauge magnets, in optical lattices. We analyze the properties of these models that make them suitable 
for quantum simulations. Within this class, we study in detail the phases of a U(l)-invariant lattice gauge theory 
in 2+1 dimensions, originally proposed by P. Orland. By using exact diagonalization, we extract the low-energy 
states for small lattices, up to 4 x 4. We confirm that the model has two phases, with the confined entangled one 
characterized by strings wrapping around the whole lattice. We explain how to study larger lattices by using 
either tensor network techniques or digital quantum simulations with Rydberg atoms loaded in optical lattices, 
where we discuss in detail a protocol for the preparation of the ground-state. We propose two key experimental 
tests that can be used as smoking gun of the proper implementation of a gauge theory in optical lattices. These 
tests consist in verifying the absence of spontaneous (gauge) symmetry breaking of the ground-state and the 
presence of charge confinement. We also comment on the relation between standard compact U(l) lattice gauge 
theory and the model considered in this paper. 
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I. INTRODUCTION 



Studying the strongly correlated regime of quantum many-body systems (QMBS) remains a very challenging task. We do 
believe that most of the interesting physics associated to them can be captured by simple models described through Hamiltonians 
involving interactions among only few of their constituents. However, even for such simple models, we are not capable of 
deriving their properties, as we still miss computational tools that are generically applicable. 

In few fortunate cases, the models constructed are exactly solvable, and almost all physical quantities can be computed exactly. 
When this is not the case, one needs to appeal to numerical simulations. A major obstacle to this approach is provided by the 
fact that the complexity of simulations increases exponentially with the number of the QMBS constituents. This means that, 
typically, one can only deal exactly with few (say up to forty) constituents. 

Some larger QMBS can still be simulated approximately. This occurs, for instance, in weakly interacting systems, where 
perturbation theory (spin-wave and its modifications), mean-field approaches and density functional theory are particularly 
successful. In all these cases, indeed, there is a small parameter - namely the interaction among constituents- that can at the 
beginning be neglected and then systematically included in the model [1,2]. 

The strongly interacting regime of QMBS is much more complex to address for the lack of such small parameter. However, 
the equilibrium properties of bosonic non-frustrated systems can be extracted using Quantum Monte-Carlo (QMC) techniques 
[3, 4]. Frustrated and/or fermionic systems are more difficult to study, but great effort has also been dedicated to them. At 
present, however, different methods often provide different results so that physical properties of these models are in most cases 
a subject of controversy (see i.e. [3] vs. [5]). 

For weakly entangled states, tensor networks algorithms (TN) are emerging as an alternative to other techniques such as QMC, 
and in some cases can also be merged with the latter to provide new numerical tools [6-10]. TN aim at providing a variational 
ansatz to describe low-energy states of QMBS, with resources that scale at most polynomially with the number of constituents. 
These techniques also give access to short-time dynamics of QMBS [11, 12], but for the moment still fail in describing their 
long-time, out of equilibrium, dynamics. 

In this situation, analog quantum simulations (AQS) appear as an alternative tool to study QMBS [13, 14]. These are experi- 
mental setups with quantum constituents (photons, ions, atoms) that are used to encode the desired QMBS [15]. State of the art 
control of these quantum systems allows also to artificially induce the desired many-body Hamiltonian acting on the constituents. 

Gauge theories are special types of QMBS, where states of the system are forced to obey local symmetry constraints. They 
stand at the core of our understanding of known fundamental interactions (electro-weak, strong, and, in a broader sense, gravita- 
tional). They are also candidates to describe the low-energy physics of frustrated antiferromagnets [16, 17]. As for other QMBS, 
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their strong coupling regime is not completely understood: in particular, they present puzzling phases, such as, confined phases 
and topological phases, which escape the standard Landau paradigm of phase classification. 

For gauge theories, and their version defined on (spatial) lattices (LGT), in the absence of fermions - or when the number of 
fermions is such that there is no sign problem- one can perform Monte-Carlo simulations. Recently, ideas about how to perform 
tensor network simulations of such theories, have been push forward as well [18]. 

Applying AQS to study gauge theories is particularly challenging since, in general, their dynamics involve Hamiltonian with 
couplings among more than the nearest-neighbors, and more than just two-body terms. This is the reason why, at present, the 
community has mainly focused on simulating the effects of static non-trivial background gauge fields on the phases of the matter. 
In this limit, the many-body interactions can be neglected, and the effect of the gauge field reduces to a specific modifications of 
the kinetic term of the matter field - hopping in the lattice- that only involves two-body interactions. 

In this context, most of the studies deal with ultracold atoms [15], and some of the theoretical proposals have been already 
realized in the experiments on the lattice [19], or, off the lattice [20, 21]. In optical lattice simulations, one can tailor the hopping 
phase/matrices properly, mainly by optical means (but also by lattice time-modulation or rotation [22, 23]). The interest in 
such studies ranges from Quantum Hall physics, where strong magnetic fluxes - very hard to obtain in condensed matter with 
real magnetic field- are needed [24, 25], to the simulation of relativistic matter and topological insulators, or extra-dimensions 
[26, 27]. 

A natural step further towards the full simulation of gauge theories with AQS, is to include the dynamics of the gauge fields. 
Such aim, for instance in two dimensions, requires the engineering of at least three-body couplings in triangular lattices and 
four-body couplings in the more common square lattices. These kind of interactions are very difficult to induce in AQS, such as, 
ultracold optical lattices [28]). There have already been some proposals to overcome these difficulties for Abelian gauge fields 
using molecular states [29], and, more recently, BECs and single atoms with internal states in [30] and [31], respectively. 

Here, we pursue an alternative solution. We propose to simulate the dynamics of gauge fields by implementing digital quantum 
simulations with Rydberg atoms (DQRS). Rydberg atoms are neutral atoms that can be excited to states close to the continuum 
spectrum, having very strong dipolar moments. Such dipoles induce long-range interactions that make possible the simultaneous 
interaction of several atoms together. The basic ideas for DQRS are described in [32] and involve the presence of ancilla atoms, 
apart from the ones entering the Hamiltonian of the system to be simulated. These are called "control" Rydberg atoms. The 
other atoms, whose interaction should encode the Hamiltonian to be simulated, are generally called "ensemble" Rydberg atoms. 
They have to be physically arranged following the pattern of the many-body interactions appearing in the desired Hamiltonian. 
For example, if one is trying to encode four-body interaction of the four atoms around an elementary plaquette of a square lattice, 
such atoms should be inside the blockade radius of a given control atom, whose dipole moment is used to implement the wanted 
plaquette interaction among them. 

One can carry out simultaneous operations on all the atoms inside the blockade radius of a control atom, operations determined 
by the state of the control atom itself. This can be achieved through a laser setup involving two- and three-photons transitions. 
In particular, one can engineer 2D lattices where the ensemble atoms are arranged on the links, while the control atoms are at 
each site and at the center of each plaquette; at same time, one can tune the lattice spacing so that all the links belonging to a 
plaquette and entering a site are simultaneously contained inside the blockade radius of the respective plaquette and site control 
qubits. With this geometry one can perform LGT simulations [33]. 

In this setup, one can perform arbitrary time-evolution with the desired (Abelian) LGT Hamiltonian, and, in some simple 
cases, one can also accomplish dissipative quantum simulations to prepare a desired state. This is, for instance, the technique 
proposed for the preparation of the ground-state of the Toric Code in [33], whose low-energy physics, in the appropriate limit, 
can be mapped to the Z2 LGT. 

In this paper, we provide a specific proposal to perform a DQRS of a U(l) LGT in two dimensions. An important ingredient 
of the proposal is the choice of a specific U(l) LGT where the constituents are qubits. Such qubits can be represented in terms 
of atoms with two logical level states, such as the ones available in standard Rydberg setups. The corresponding formulation of 
U(l) LGT is also known as link model or gauge magnet [34-37]. 

The outcomes of our proposal are multifold. First, the alternative formulation that we consider is not only more practical 
to be implemented - the local Hilbert space is finite- but also - at least it is our belief- closer to the language of the QMBS 
community than the original Hamiltonian formalism proposed by Kogut and Susskind [45]. Within this general framework, we 
show how to derive the specific U(l) gauge magnet we are interested in. 

Furthemore, we introduce an exactly solvable version of gauge magnet, where we are able to provide the exact ground-state 
wave function. We use this example to illustrate two paradigmatic properties of gauge theories, i.e., the absence of spontaneous 
breaking of the gauge symmetry, and the presence of charge confinement. The quantitative check that such properties hold can 
be used by experimentalists as a sine qua non test of the quality of the simulator. 

Finally, we determine a specific protocol to prepare the ground-state of the model and perform arbitrary, out of equilibrium, 
dynamical simulations. As well known, one can extract the excitation spectrum of the model by performing these simulations. 

Our approach is relevant because it allows the simulation of a LGT with continuous symmetry groups away from the special 
points, as studied in [33]. It provides a generalization of the techniques described in [38-40], such to permit to prepare ground- 
states of Hamiltonians that are not frustration-free. Furthermore, we explain how to perform out of equilibrium long-time- 
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evolution of the U(l) gauge magnet with a DQRS, following the prescription of [33]. 

It is worth to notice that our ground-state preparation scheme successfully merges several available techniques. It is, indeed, 
made of two steps. First, by taking advantage of the fact that the square lattice is bipartite, we decompose the Hamiltonian of 
the gauge magnet as the sum of two Hamiltonians, each of them for separated is frustration-free. This allows us to prepare the 
ground-state any of the two with standard dissipative protocols. The second step involves a driven evolution via DQRS to the 
final Hamiltonian of gauge magnets. The turning on of the part of the Hamiltonian that has been neglected in the first instance, 
can be adiabatic or not. Indeed, the efficency of the ground-state driving can be improved by following the chopped random 
basis (CRAB) protocol [41, 42], which allows for closed-loop optimization experiment, i.e., the preparation time is optimized 
recursively by running the DQRS. In this way, we manage to prepare the ground-state of the complete Hamiltonian, a very 
complex, strongly interacting, non frustration-free, i.e., frustrated Hamiltonian. In practise, such preparation can be applied to 
any system governed by a local Hamiltonian in bipartite lattices. 

The paper is organized as follows. In section II, we introduce the concepts necessary to formulate a LGT In doing so, we 
assume that the reader is familiar with many-body quantum systems, while the knowledge of continuous or Lagrangian formula- 
tion of gauge theory is not assumed. In particular, we perform a constructive derivation of the Hilbert space and Hamiltonian of 
an arbitrary Abelian LGT (this section is complemented by the examples contained in the Appendix). In section II B, we review 
the properties of the U(l) gauge magnets as a specific instance of generic Abelian LGT using the language and the notation 
introduced in the previous section. We discuss how to split the Hamiltonian in the sum of two frustration-free Hamiltonians. 
This is a crucial step both for applying the QDRS protocol of the ground-state preparation, and for analyzing the absence of 
spontaneous symmetry breaking and the presence of charge confinement. We propose these phenomena as key experimental 
observables to probe whether the simulator is working properly. In section III, we start by reviewing the ideas beyond digital 
quantum simulation and the mesoscopic Rydberg gate [32]. We discuss in detail the setup necessary to perform a DQRS of the 
U(l) gauge magnet time-evolution, which permits the study of long-time, out of equilibrium, dynamics. We conclude by giving 
special emphasis on how to prepare the ground-state of the U(l) gauge magnet with a mixture of dissipative and adiabatic evolu- 
tion and discuss the relation between our model and standard compact U(l) LGT. We propose two possible sets of measurements 
that can be used by the experimentalists to check both absence of spontaneous symmetry breaking and charge confinement. For 
reader's convenience, the main features of the experimental set-up and the experimental sequence are summarized in section 
HIF. 



II. CONSTRUCTING AN ABELIAN LATTICE GAUGE THEORY 

We introduce an Hamiltonian formulation of Abelian lattice gauge theories. Gauge theories are generally studied in the 
Lagrangian formulation, by introducing an action invariant under local transformations, and quantizing it via the path integral 
formulation of Quantum Mechanics [43, 44]. By discretizing the space-time, one obtains a lattice description that can be treated 
with the tools of statistical field theories. In this formulation, the corresponding quantum Hamiltonian can be derived from the 
transfer matrix of the model, as as originally proposed by Kogut and Susskind [45] (in the same way one can map the classical 
2D Ising model at finite temperature on to the quantum Ising chain at the zero temperature with a transverse field [46]). 

Here, we take an alternative route. Indeed, an Abelian lattice gauge theory is a particular case of a many -body quantum system, 
where the constituents are arranged on a lattice and the states (and the observables) are invariant (as we will explain below) under 
local transformations, which are elements of a given Abelian group G. Starting by this operative definition, we can construct an 
Hamiltonian for a LGT. It is worth to notice that, for a given G, infinitely many gauge invariant theories - and Hamiltonian- can 
be defined, although many of them may be equivalent under renormalization group flow, i.e., may share the same fixed point and 
critical properties. Baring this consideration in mind, the basics steps to follow to construct a gauge invariant Hamiltonian for 
a lattice system, are: i) the choice of the gauge group G that, as we will see, has immediate consequences on the definition of 
the Hilbert space of the model; ii) the definition of the operators defining the gauge transformations; and iii) the determination 
of the set of operators that can be used to build the Hamiltonian, which has to be invariant under gauge transformation. This 
constructive approach should sound more familiar to people working in the field of quantum many-body systems and condensed 
matter rather than the one originally proposed by Kogut and Susskind, as it is not requiring a prior knowledge of gauge theories. 
It also makes explicit that each of the points i)-iii) entails a large freedom. The original Kogut-Susskind proposal reveals itself 
as a specific choice of LGT. Also, using this approach one can interpret different LGT known in the literature as link models 
or gauge magnets [34-37] as a truncated version of the original Kogut-Susskind proposal. We would like to mention that this 
approach is rooted in the celebrated works by Kitaev about the Toric Code and Quantum Doubles [47] that we follow closely. 
When defining the steps i)— iii) we will explicitly construct the simplest Z2 LGT, so that the reader will have an immediate 
example on which she/he can understand the more general abstract formulation. Several other examples are contained in the 
Appendix. 
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Figure 1: a) A lattice gauge theory is many-body system defined on a lattice jSf. The constituents |/?)/ are attached to the links Z. The lattice is 
oriented, meaning that links inherit the standard orientation of the embedding space, depicted in the figure as arrows. When defining operators 
acting on the states attached to links, one has also to specify an orientation, b) and c) The building blocks of a LGT include There two sets of 
four-body operators, b) The one attached to sites s, A s , which we also refer to as star operators, are used to impose the gauge invariance. c) 
The ones attached to plaquettes p B p , which we also refer to as plaquette operators, describe the dynamics. We define both A s and B p by using 
a customary oriented star and plaquette. 



A. Gauge invariance, the gauge invariant Hilbert space and operators 



i) The Hilbert space. In order to construct a LGT one first need to specify the symmetry group G. For example, the simplest 
choice is to consider G = Z2, which is the group formed by two elements, l,e, with the following multiplication table 1 X e = 

ex\ = e.e x e = I x \ = I. 

The constituents of the LGT are attached to the links / connecting sites s of an oriented lattice Jzf. The fact that the lattice is 
oriented means that, when addressing a link, one has not only to specify its position but also its orientation (the direction to follow 
in order to walk along it, i.e., in 2D either from left to right or from up to down or viceversa). The need for considering oriented 
lattices can be understood through the physical picture of the links as currents moving from a site to another, the orientation is 
the direction of the current (and the gauge condition is the analogue of the conservation of the current at each site). 

Each of the constituents attached to the oriented links is described by a vector in a local Hilbert space fj. If the lattice has N 
links, the global Hilbert space is defined as the tensor product M' = Y,® N ' . 

The choice of ~fi is not unique, and this is the first freedom at hand in the construction of LGT. In particular, the original 
formulation by Kogut and Susskind uses as local Hilbert space the group algebra C(G), so that the space of constituents is the 
vector space generated by linear combinations of elements of G with complex coefficients. To each element of the group G, 
one associates a vector \g), with the property (h\g) = 5„/,. Here 5 is the Kronecker delta. In the simplest case of Z2, the 
local Hilbert space is generated by the two orthogonal vectors |1) and \e), and it is isomorphic to C 2 . In the Kogut-Susskind 
formulation, the dimension of the local Hilbert space is strictly related to the number of elements in the group. Hence, it becomes 
infinite for continuous groups like U(l). However, as we will show below, it is perfectly consistent to consider a local Hilbert 
space whose dimension does not depend on the number of elements G. 

As an example, the space of constituents of a 2D lattice is sketched in Fig. 1. Since we are dealing with an oriented lattice, 
we also orient the states. Following physical intuition of a links as currents between sites, the change of orientation corresponds 
to an inversion of the currents. With the standard choice of the local Hilbert space as group algebra, "fi = C(G), the physical 
intuition translates in associating the inverse in G to a flip of the orientation. This can always done since the states are elements 
of a group that by definition have their own inverse. Hence, changing the orientation of a link is equivalent to the operation 

-H/T 1 }-/. (1) 

The detailed graphical explanation of the role of the orientation is given in the Appendix B, Fig. 1 1 iv). It is worth to observe that, 
for the guiding example of G — Z2, the orientation plays no role since e _1 = e and 1 = 1, hence, \e)i — \e)-i and 1 1)/ = 

The only requirement on the local Hilbert space, needed to construct a consistent LGT, is that fi has to be isomorphic to a 
representation of the symmetry group G. In particular, the choice of C(G) induces the choice of the fundamental representation 
of the group, Rf(G). The action of an element of X(g) g Rf(G) on \h), is defined as 

X(g)\h) = \(gh)), (2) 



This means that the matrix representation of X(g) reads 

X(g) 



= E \gh)(h\. 

heG 



(3) 
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For G = Z 2 , it follows thatX(l) = |1)(1| + \e)(e\ andX(e) = \e)(l\ + |l)(e|, which coincide with the identity matrix and the o x 
Pauli matrix, respectively. 

When we change the orientation of a link, accordingly to (1), we also have to invert the operator that is acting on it, 

\{gh))i \{h- l g~ l )U =X_,(A" l g- l h)\h- l )-i =X_ l (gf\h- l )_, (4) 

where / is an arbitrary oriented link of Jz? and we have used the fact that, for Abelian groups, h~ 1 g~ 1 h — g , and, for unitary 
representations, X(g~ l ) = X(g) '. The irrelevance of the orientation in the Z 2 case translates in that both X{\) and X(e) are 
Hermitian. 

ii) Gauge invariance and the physical Hilbert space. A generic gauge transformation is defined by choosing a group element 
for each site of the lattice and rotating all the links entering that particular site with the matrices representing the rotation for 
that particular element. The building blocks of gauge transformations are thus obtained by acting on all links entering a given 
site with X(g). The operator that induces such rotation is often referred as star operator and for the specific case of a 2D square 
lattice is defined as 

A s(g)= x i(s)-\hh,hi 2 ,hi 3 ,hi 4 )^\gh h) gh h ,ghi 3 ,gh k ), (5) 
le{hh 

where s is a generic site of the lattice and {/,}, i = 1, ... ,4, is the set of all the links entering s, as depicted on Fig. 1 b). In view 
of (1) and (4), this operator acting on the links of a lattice oriented in the conventional way (such as in Fig. la) becomes 

A s(g) ■ \h-i^h-i z ,h l3 ,h k ) -> \ h i { 1 S~ l ^2 g~ l i8 h h'S h h) ' ( 6 ) 

where the links are numbered clockwise starting by the one on the top. In terms of X(g) operators, the above expression is 
equivalent to A s (g) — Xj| (g)xf (g)X/ 3 (g)Xi 4 (g) . The choice of the group element g can vary from site to site so it would be more 
appropriate to use g s but in order to simplify the notation we stick to g. For G = 1%, the only non-trivial A s operators are those 
associated with the group element e and are defined as A s (e) = (7^ (7^(7^(7^. The action of A s {e) then induces spin-flips on all 
qubits entering the specific link s. 

States are invariant under the gauge transformations at site s if they are eigenvectors of the A s (g) operators with eigenvalue 
+ 1, 

A s (g)\ ¥ ) = \ ¥ ),Vg€G. (7) 
Since generic gauge transformations are product of local ones, a state is gauge invariant if 

&({g})\Y) = ®Mg)\v) = \Y),VgeG. (8) 

This leads to the definition of the physical Hilbert space, M'g, as the set of those states in H that are gauge invariant, 

Jtb = {\Y)} s.t. n{g})\Y) = W),VgeG. (9) 
It is important to notice that the physical Hilbert space is a subspace of the original Hilbert space, 

3tf G e = r m ~ &W® N , (10) 

For G = Z2, the requirement of invariance at a specific site under the action of A s (e) reduces the dimension of the Hilbert space 
from 16 to 8. In order to determine the allowed states, it is better to diagonalize o x , G x = |+)(+| — |— )(— |. In this basis, all those 
states formed as tensor product of four eigenstates of o x with an even number of | — ) are eigenstates of A s (e) with eigenvalue 
+ 1, and thus gauge-invariant. As we will describe in the following, one can give a nice geometrical interpretation to those states 
in terms of closed string of |— ). 

iii) Operators compatible with the requirement of gauge invariance. Once defined the physical Hilbert space, we focus on the 
operators that are compatible with the requirement of gauge invariance. By definition, the LGT Hamiltonian H has to respect 
the local symmetry, i.e., commute with all star operators A s (g s ), 

[H,A S ] = 0, VA S . (11) 

Our goal is to built up H as a sum of local terms. 

To this aim, we introduce a set of operators Z(r), acting on the tensor product of ij^Yg. Again, there is a lot of freedom 
in defining the factor y g , since the only requirement is to support an irreducible representation (irrep) of the gauge group (for 
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a definition of irreducible representation refer to [48]). When the local Hilbert space is the group algebra, the operators Z(r) 
acquire the following form 

Z(r)='£R r (h)®\h)(h\, (12) 
heG 

where R r (h) is the matrix representing h in the irreducible representation r of G. Since all irreducible representation of Abelian 
groups are one dimensional and isomorphic, for Abelian theories we can drop the index r from R r and think of them as acting 
only on "fi since, in this case, R(h) is just a phase. For example, in the case of G = 1%, the only non-trivial choice for Z is 
Z = o^ = \l)(l\-\e)(e\. 

The relation between Z and X(g) is encoded in the commutation relation 

ZX(g)\h) =R{gh)\gh) =R(g)XZ\h) (13) 

This immediately suggests a minimal choice to fulfill (11). We can, indeed, consider as building blocks for H the product of 
four Z operators acting on links around elementary plaquettes of the lattice 

B p = Z,(r), (14) 
le{k} P 

that we generally call plaquette operator . Here, {li} p , i= 1 , ... ,4, is a set of links belonging to the plaquette p and anti-clockwise 
oriented starting from bottom, as sketched in Fig. lc. By rewriting B p for the standard-oriented lattice, as in Fig. la, we obtain 

B p =Z h (r)Z h (r)Z(r)\z(r)\. (15) 

In the elementary case of G = 1%, the plaquette operator is just B p = of of of of . 

As a plaquette and a star operator share none or two links, / and it is sufficient to verify the relation 

[Z l (r)®Z l ,{r),X l {g)®X}(g)]=Q, (16) 

in order to check that B p and A s commute. As the above relation is a direct consequence of (13) (and that for any representation 
R{g~ l ) = R~ 1 {g)), the desired result 

[B p ,A s (g)]=0, V{p,s}GJ?,g€G, (17) 

holds. 

As obvious for an Abelian gauge group - the case of interest in this work- any operator X[(h) commutes with A s (g s ). Hence, 
it follows that any hermitian functional of B p and Xi(h) is a good gauge invariant Hamiltonian for Abelian gauge theories. In 
particular, we focus on the linear combination 

H(6) = -cos9Y,B p + smeY,Xi{g)+H.c, (18) 

p I 

where p are the elementary plaquettes of the lattice and / are the links. This coincides with the Kogut-Susskind Hamiltonians 
[45]. In the specific case of G = Z2, it can be written explicitly in terms of Pauli matrices H(6) = — cos 9 Y. p Yhep °f + 
sindZi of {g). 

On physical ground, the B p and X[ operators describe magnetic and electric interactions, respectively. For the specific example 
of G = Z2, both electric and magnetic fluxes can only have two values, ±1. 

It is worth noticing that the presence in the Hamiltonian of terms involving four-body interactions is a direct consequence of 
gauge invariance. There have been several attempts to obtain a gauge invariant Hamiltonian starting from a model with only 
first neighbour interactions. For instance, by defining two-body Hamiltonians on a coarse grained lattice [49], one can obtain 
in some case equivalent Hamiltonians to those with four-body interactions. However, at the moment these type of constructions 
are restricted to exactly solvable models, and it is unclear how to generalize them to a generic Hamiltonian, as the one we will 
consider in the following. Thus, in this paper we stick to the idea that gauge invariance (on square lattices) requires four-body 
interactions. 



B. Gauge magnets or link models 



Up to this point, we have seen the general requirements for obtaining a LGT That is to say, a LGT is characterized by i) the 
presence of a gauge symmetry that permits to ii) identify a subspace of the full Hilbert space as the physical Hilbert space where 
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iii) only gauge invariant operators can be considered. In order to get the original Kogut-Susskind proposal, we need to make 
specific choices for all the three points. Indeed, for i) we use as a local Hilbert space the group algebra C(G). This immediately 
implies that in ii) the X(g) operators are obtained by considering the regular representation of the rotation by an element g of 
group G. Furthermore, in iii) we have constructed the Hamiltonian using the simplest closed path of the lattice, the plaquettes. 
Nevertheless, any of the above choices may be changed, resulting in a different gauge theory. Here, we focus on alternative 
choices to i). The reason is that the C(G) algebra for a continuous group is infinite dimensional, while the implementation we 
have in mind of LGT simulator with optical lattices can only deal with finite dimensional local Hilbert spaces (see section III). 

Our aim is to construct an Abelian LGT with the smallest possible local Hilbert space, independently on the number of 
elements of the group G. This leads to LGT that have been called gauge magnets or link models [34—37] in the literature, whose 
particular case is that U(l) LGT with finite dimensional local Hilbert space. 

Let us start by comparing those LGT to the formulation of Kogut-Susskind. In order to do so, let us consider again the 
(possibly infinite dimensional) local Hilbert space C(G). A celebrated theorem of group theory [48], allows us to truncate C(G) 
to a finite dimensional Hilbert space, without destroying the gauge symmetric structure. The theorem states that the regular 
representation (the one that acts on C(G)) can be written as the direct sum of all possible irreducible representations (with a 
multiplicity equal to their dimensionality). 

As a consequence of this theorem, we are guaranteed that there is a change of basis such that the operators X(g) become block 
diagonal, 

X(g) = ® r j:R r (gy- i \(r,i))((rJ)\, (19) 

ij 

where we have explicitly written the block structure of X(g) labeled by the r irrep, and the R r (g) u] are matrix elements of the 
irrep r of g € G, i,j — 1 • • -dim(r). We call the rotation matrix that brings all X(g) to the block diagonal a((r,i),g), so that we 
can express the Z(r) operators in the new basis as 

Z(r)= £ z(r)WM\{p,i))({q,j)\. (20) 
{p-')M,j) 

where the matrix z(rp p '^'^' is z{r)^'^> — T 1 g a ((Pi')ig)Rr(g)<^ l (g, We call the new basis with abuse of notation 

{\r)} basis. 

In the {\r)} basis, we can safely truncate the local Hilbert space of the LGT without any effect on the symmetry requirements. 
Indeed, we just have to include at least one of the diagonal blocks of X(g). However, it could happen that by keeping just one 
block either X(g) or Z(r) become trivial (that is indeed the case for Abelian theories). For this reason, we need to keep at least 

two irreps. In this case, the local Hilbert space will have dimension d = dim(r\) H dim(r n ), where one can stop at the first n 

that provides both non-trivial symmetry requirement, and non-trivial dynamics. 

In the following, we will provide the specific example of the t/(l) gauge magnets, since it is the one we are interested in 
simulating with optical lattices. In the Appendices, we provide further examples for generic Zn gauge magnets. 



C. U(l) Gauge magnet 

Using the ideas of the previous section, we can easily study the U(l) gauge magnet. A gauge magnet differs from the standard 
LGT by the choice of the local Hilbert space. Rather than using the group algebra, one retains at least two arbitrary irreps, 
r\ (Br2- Since U(l) is Abelian, its irreps are one-dimensional: it means that we can work with qubits. We could choose any 
irreps r\ and ri, the simplest choice is to select the trivial irrep where all the element are mapped to 1, and the irrep where a 
given element g SU(1), distinct than the identity 1, is mapped to the phase e ia s. As we are interested in faithful irreps, a g has to 
be chosen not commensurable with 2n, i.e, ^ {Rationals}, if g is s.t. g" ^ 1, Vn € N. Note that asking local invariant under 
a g as such it is sufficient to ensure invariance under any local U(l ^transformation. We denote the eigenvectors of the X(g)s as 
|+) and |— ). This means that the X operators are 

X(g) = \+)(+\ + \-)(-\e ia s, (21) 

while the Z operator reads 

Z=\+)(-\. (22) 

This specific model can be mapped (on bipartite lattices) by a unitary transformation to the one studied in [50]. 

In order to describe the physics of the model, it is important to notice that we can characterize J^q graphically. In the |+) and 
| — ) basis and with the standard 2D orientation of Fig. 1, the gauge condition (7) on a given site s selects among the 2 4 states of 
the four links entering s only six possible configurations 



|+, +,+,+>, I- +,-,+),!+, -,+,-),!-,-,- -),!-,+, +,-),!+,-,-,+). 



(23) 
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These are shown in Fig. 2. Similarly to standard LGT, we can interpret the U(l)-invariant Hilbert space as the space of closed 
strings formed by |— ) states, onto the vacuum of |+) states. Indeed, as any of the six allowed states contains an even number of 
|— ), at each site of the lattice, for any ingoing |— ) there is also the corresponding outgoing one. Hence, the string cannot end 
on a (bulk) site, i.e. only closed string (or string touching the boundary) are compatible with gauge invariance. In addition, it is 




b) 



Figure 2: a) The requirement of gauge invariance of the states using (7) selects only 6 out of the total 16 states of the four qubits entering a 
given site. In the standard 2D orientation, one can visualize the 6 allowed states by coloring in red states |+) and blue |— ). b) If we consider 
the standard Z2 LGT the gauge symmetry condition of (7) selects 8 out of the 16 states of the four qubits attached to the links entering a given 
site. They include the same 6 states than those of the U(l) gauge magnet of panel a). The two extra states are represented here, and are those 
responsible of the existence short closed loops. 



worth to notice that the six allowed configurations do not allow to construct a closed string of finite size (without reaching the 
boundary). This means that, for lattices with periodic boundary conditions, strings can only close by wrapping around the whole 
lattice. This implies that their length is at least equal to the lattice size. This is an important difference with respect to standard 
LGT. In standard LGT, indeed, closed strings can be of arbitrary length, the shortest being the strings around a single plaquette. 

In order to clarify the origin of this discrepancy, one can consider the Kogut-Susskind Z2 LGT introduced in the previous 
section. There, the group algebra is two-dimensional, that is, the local Hilbert space is still made of qubits. As we have already 
shown, the gauge condition selects 8 out of the 16 states of the four links entering a given site. It turns out that 6 of them coincide 
with the one of the U(l) gauge magnet of Fig. 2 a), but there are two extra states 

I- -+,+),]+,+,-, -), (24) 

as shown in Fig. 2 b). 

These states are, indeed, the ones needed to close a small loop of |— ) (for example a single plaquette). In fact, this small 
difference produces completely different strings patterns, as we can appreciate in Fig. 3. 

The different string pattern have deep consequences on the physics of the models. Standard LGT, indeed, have at least two 
phases, a confined phase and a deconfined phase. In the confined phase, short closed string abound, while long closed string are 
very rare and vice-versa in the deconfined phase [51]. In the 1/(1) gauge magnet, by just noticing the absence of short closed 
strings, we already have a strong indication that the phase diagram of the model is very different from the one of the standard 
t/(l)LGT. 

Let us analyze it in details. First, let us consider the dynamics induced by the Hamiltonian (18) when 0=0. Contrary to 
what happens in standard Abelian LGT, for the gauge magnets two plaquettes operator sharing one link are not commuting, i.e., 
generally [B p ,B' p ] 7^ 0, . This means that the ground-state is not simultaneous eigenvector of all the Bp's. However, for the sake 
of the present discussion we can, as suggested in [50], circumvent this problem on bipartite lattices by considering, first, only 
half of Hamiltonian, and, afterwards, the effects of the other terms in (18). We start with 

flb = (25) 

Py 

where with p y we label half of the plaquettes, the ones drawn in yellow in Fig. 4. In this way, Ho only contains operators that do 
not share any link, and, thus, all commuting. 

The model (25) is exactly solvable, and we can write its ground-state as the action of a series of projectors onto a given 
reference state. It is important to notice that the reference state should be i) gauge invariant, and ii) not belonging to the kernel of 
the various B p considered. Indeed, in contrast to what happens for standard LGT, the B p here has only two eigenvalues different 
from zero, equal to ±1. It is easy to check that a possible reference state fulfilling these requirements is the one depicted in the 
central panel of Fig. 4, consisting of two strings each of length L 2 /2 wrapping around the lattice of size Lx L with periodic 
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Figure 3: Here, we depict a 4x4 lattice with periodic boundary conditions in both directions. Red qubits are in the state |+), while blue qubits 
are in the state |— ). The blue qubits form closed string pattern as required by gauge invariance. On the left, we show a possible gauge invariant 
state of the standard Z2 LGT. On the right, the same plot is presented for the U(l) gauge magnet. In both cases, the states can be mapped to 
closed string configurations. There is an important difference between the two cases, however. In Z2 LGT, closed strings can form arbitrary 
small closed loops. In the U(l) gauge magnet, due to the absence of the states in Fig. 2 b), strings are forced to close by wrapping around the 
whole lattice. This means that their typical length exceeds that of the linear lattice size. Blue lines are drawn as a guide to the eye to recognize 
the closed string wrapping around the periodic boundaries of the lattice. 
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Figure 4: a) Lattices with an even number of sites and periodic boundary conditions are bipartite. This allow to consider the model with half 
of the plaquettes (filled in yellow in the figure) turned on. In this way, the different yellow plaquettes do not share any link and the U(l) gauge 
magnet Hamiltonian is built up of mutually commuting terms. This implies that the ground-state is simultaneous eigenvector of all the terms 
appearing in the Hamiltonian. b) Reference state used to construct the ground-state of the Hamiltonian (25) through the action of projectors 
(26). c) An example of reference state used to obtain an excitation with gap A = 2. 



boundary conditions. We call this reference state | yb) • We can define the three different projectors on to the subspace of different 
eigenvectors of (B p +B* p ) as 

P* = ±(B p +Bl)(B p +Bl + l), (26) 

P; l = \(B P +Bl)(B p +Bl-t), (27) 

P° p = - (B p + B\ + 1) (B p + B\ - 1) , (28) 
so that the ground-state |£2(0)) is proportional to 

l^(0))-®<IVo). (29) 

Py 

By construction the above state is not null and minimizes the energy as 

1^(0)) - |l) Pv , (Bp, +<) |l) Pv = \l) Py ,V Py G {yellow plaquettes}. (30) 

Py 
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As, in principle, all the eigenvectors of Hq can be constructed by applying the projectors (26-28) on an appropriate gauge 
invariant state, it follows that there is a gap to the first excited state. Such gap is at least equal to one, corresponding to a state 



m{o,p'})) = 



|l> ft ®|0). 



Py^Py 



B p'y+ B p>, 



|0U=0. 



(31) 



However, it is easy to argue that such a state cannot satisfy the gauge invariance requirement on each of the sites at the vertices 
of the plaquette p' . Hence, |O(0)} describes a gapped phase with gap 2. The simplest excited state with such energy is 



mi-hp'})) 



Py^Py 



/Py> 



B, 



-B 



(32) 



as \<i>({-l,p' y })) oc <g) Py ^ Py P py |Vo)- ^ is worth to notice that first excited states of the form \(S>({0,p' y }, {0, p"})) = 

®p ^p' p" \ ® \®)p' ® \®)p" a ^ so ex ist and are gauge invariant, but cannot be obtained just by applying projectors on the 



reference state \ \jf ) as P2 \ Yo) 



. ,, ril , — 0, \fp y in the yellow plaquettes' sublattice. In this case, the procedure is more complex and it 
involves a B operator acting on a plaquette of the other sublattice, as illustrated in the right pannel of Fig. 4. It is found that 



|<&({0,/ },{0, Mo }))-«Bo|n(0)) 



(33) 



where do, ro, mo, h are the plaquettes in the yellow sublattice that surrounds the plaquette of the complementary sublattice. 

In this simple model, we can compute exactly how the presence of static external charges modifies the ground-state. A static 
charge ±1 at site s modifies the gauge condition (7) to 



A s (g)\v) =sxp(±ia g )\Y). 



(34) 



That is to say, the presence of static charges modifies the allowed string configurations. Open strings are, indeed, allowed if they 
start and end on the charges. This, clearly, affects the properties of the ground-state of the system. With two charges ±1, the 
ground-state is orthogonal to the one without static charges. This implies that we have to change both the reference state and the 
set of projectors, in order to construct the new ground-state 



|fl(0) d 



)P Py \Y±i) 



(35) 



A simple candidate for |v±i) can be obtained by transforming one of the two closed strings of |— ) contained in \y/o) into an 
open one. As illustrated in in Fig. 5 for a 4 x 4 lattice with periodic boundary condition, this can be done by flipping one or 
more consecutive |— ) (blue) links to (red) |+) links. This creates two static charges of opposite charge at the two ends of the 
blue string, where modified gauge condition (34) holds. Two possible choices of the reference state y±\, which differ on the 
position of the static charges, are sketched in Fig. 5. The +1 static charges are denoted by filled black dots, while —1 ones by 
empty dots. 

Now, the projectors P p are determined by the requirement of minimizing the energy while respecting the new gauge condition, 
which depends only on the position of the ± 1 charges. Since the above defined | y/±i) differ from | Yo) on ly by f ew flipped links, 
located between the two ends of the open string, the projectors P Py may be distinct from the Pt only for the plaquettes p y (of the 
yellow lattice) interested by such flips, p y G {flipped region} (pale yellow shaded plaquettes of Fig. 5). In fact, as manifestation 
of orthogonality of the ground-state with back-ground charges to the ground-state without charges, |y/±i) is annihilated by the 
projectors (26) acting on the plaquettes between the two ends of the open string. Hence, the proper choice of projectors in such a 
region, which minimizes the energy and does not annihilate the state, is given by the projector of (28), P Py = P2, Vp y & {flipped 
region}. 

As a consequence, the ground-state energy of the system with two static charges is higher than the one without the two charges. 
The energy gap A q is equal to the number of P° projectors, i.e., the number of yellow plaquettes that contain flipped links. Up 
to artifacts of the discretization, such number is proportional to the number of flipped links itself, rip , i.e., is proportional to the 
distance r, in lattice spacing unit, between the charges, A q {r) oc r. As it can be easily deduced from Fig. 5, the exact relation 
between the gap and the number of flipped links is 



/if 

Aq = n F -Int[y], 



while the distance between the charge depends on rif as 



rip 



rip 



» F -Intbf] (1, )+Intbf](0,l) 



2Int 



np 



] 2 -2« f Int[^] 
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Figure 5: The presence of static external charges modifies the ground- state, a) Sketch of the reference configuration used to obtain the 
ground-state with two static first neighbors charges (presented as filled black dots for charge +1 and empty dots for charge —1). The blue 
string os states |— ) is broken in between them. The plaquettes states along the broken part of the string are orthogonal to the eigenstates of 
Bp with eigenvalue 1 . These plaquette are filled in with a pale shade of yellow in the drawing while regular plaquettes that contribute — 1 
to the ground-state energy filled in solid yellow. This implies that the ground-state energy in the presence of static charges increases by an 
amount proportional to the inter-charge separation, manifestation of charge confinement, b) Sketch of the reference state for charges at distance 
r = 2v2. The reference state induces a further pale yellow plaquette with respect to the reference state on the left and has thus higher energy. 



where Int[jc] is the integer part of x. The above formula is obtained by taking in account the "zig-zag" behavior the displacement 
of the charges shows when one additional link is flipped. 

These two relations together are the footprint of charge confinement, a phenomenon that the present model shares with many 
other gauge theories, among them, QCD. This means that the gap scales as 

\(r) = or, (36) 

where a is called string tension. In the long distance regime, i.e., np ^> 1, it follows that A c/ ~ ^ and r ~ hence, A q ~ -^=, 
and the string tension for this model is a ~ . 

At this point, we consider the effect of the link term for 0^0. By the Hamiltonian cos0//o + sin0 Y.i{Xi(g) +Xj(g)). 
The last term of the Hamiltonian, cos 9 Ho + sin0 +xj(g)) is certainly dominant close to = |. In this phase, the 

ground-state is the product state 



in 



(f)> = ®l->'. < 37 > 



regardless of the a ? we have chosen to represent a generic 1/(1) element g. Hence, one expects to encounter a phase transition 
when grows from zero to f . We observe that i) 



tfoin(f); 



= 0, (38) 
i.e., the expectation value of the plaquette part of the Hamiltonian vanishes in the ground-state at 6 = ?, and, as well, ii) 

£(Xi+X/)|l) p = 4(l+cos(%)|l)„V/> J (39) 

lep 

which implies + x/)|n(0)) = (1 + cos a s )L 2 \Q(0)). From i) and ii), it immediately follows that the states |fl(0)) and 

|n(f )) are not deformed, but simply shifted in energy, by the change of 0. Hence, the transition between them is a first 
order phase transition (level crossing) where the expectation value of any (yellow) plaquette operator can be taken as an order 
parameter which jumps abruptly from — 1 to 0. 

Furthermore, the relations (38,39) allow to compute the critical value of the coupling 9 C simply by equating the energy per 
plaquette of |Q(0)) and |n (f )) 

/ j cosf cc ) \ 

9 C : -cos0 c +4sin0 c (l +cosa g ) = 8sin(0 c )cos(a g ), -> 9 C = arctan 8J . (40) 



Also in the product state phase, static charges are confined. The minimal configuration containing to two charge excitations of 
opposite sign, ±1, is given by the shortest open string of |+) -links connecting the charges. These configurations have an energy 
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Figure 6: Scaling of the gap A between the ground-state and the first excited state of the Hamiltonian of (41) as a function of the coupling c w 
for a 4 x 4 lattice with periodic boundary conditions in all directions. The gap decreases and eventually saturates to its finite size value. On a 
infinite lattice the results of [50] predict that it would vanishes at the point at c w = 1. 

cost per link that is constant and equal to 2 — 2cos(a ? ). Again, the linear behavior of energy gap with the charge distance is a 
manifestation of the confinement of charges. 

Finally, we are ready to study the whole Hamiltonian (18). We notice that although the actual shape of ground-state in the 
plaquette dominated phase |£2(0)) is modified, the conditions (38-39) are not, as the plaquette ground-state has always an equal 
number of plus and minus links. This implies just a change in the actual value of 6 C , but not in the nature of the phase transition 
that remains a first order level-crossing phase transition. 

Hence, all the properties of the plaquette dominated phase can be studied at = 0, as a change of 6 only induces a shift in the 
energy of such state and of the states associated to plaquette excitations. 

In order to obtain the full gauge magnet Hamiltonian, we can add the other half of the plaquettes to the Hamiltonian (25) 
adiabatically 



where p w are the white plaquettes in Fig. 4. By varying c w from to 1 the system is driven to a gap-less phase, [50]. The 
dependence of the gap on c w for the 4x4 lattice is shown in Fig. 6. There we appreciate that the gap systematically decreases 
and eventually saturates to its finite size value. 

We conclude the section by making a comparison between the regular and the two-dimensional representation of U(l) LGT. 
The truncated formulation has the clear advantage of being easier to handle numerically, and to simulate in optical lattice. The 
relevant question is whether the truncated and the non-truncated theories share the same quantum properties. It is worth to notice 
that the gapless phase encountered here has no analogue in the standard compact U(l) LGT, where the only appearing phase is 
the gapped confining phase [52]. 

There are several interesting physical phenomena that are however common to both models. Here, we just list two of them, 
which we propose in the following section as possible test to check whether the desired physics is correctly realized in the quan- 
tum simulator. The first one is related to an important theorem that states that gauge symmetries can not be broken spontaneously 
[46]. The physical consequence of this fact is that only gauge invariant operators can acquire a non-vanishing expectation value. 
This implies that the expectation value of, i.e., a single Z operator on a link has to be zero, property that can be easily checked 
with current experimental resources. The second one is related to charge confinement. In the U(l) link model, we have seen that 
both phases exhibit charge confinement. This is similar to the charge confinement encountered in the Kogut-Susskind formula- 
tion of compact 1/(1) LGT. As discussed in section HIE, the presence of charge confinement can be verified experimentally. 

Definitely, in order to establish the correct relation, beyond analogies, among the standard compact QED and the present link 
model - including a full understanding of the renormalization group flow of the U(l) gauge magnet- further studies are needed. 



(41) 




D. Exact diagonalization and Tree Tensor Network of the U(l) gauge model on small lattices 



With a specific choice of open boundary conditions the model is exactly solvable. Here, we performed an exact diagonalization 
of a system with periodic boundary conditions in both directions to see which of the features of the exactly solvable model survive 
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Figure 7: In order to perform exact diagonalization on the 4x4 lattice, we need to appropriately reduce the number of degrees of freedom 
by using the constraints induced by gauge invariance. There are several way to perform the reduction. Here, we introduce an RG inspired 
procedure that can be generalized to larger lattices and would immediately provide a tree tensor network ansatz for symmetric states (see also 
[18]). The procedure consists in identifying blocks of four spins and project only on those configurations selected by the gauge conditions 
(those sketched in Fig. 2). This produces a new lattice of effective sites with dimension 6. We can iterate the procedure and obtain a final layer 
where only few spins are left. At this point, the effective Hamiltonian can be diagonalized directly. If one needs to consider larger lattices, the 
projection onto gauge invariant states has to be followed by a projection onto low-energy states, such that one can iterate a larger number of 
time the RG steps. See, e.g., [18, 53]. 



the presence of periodic boundary conditions. Clearly, with this approach we are limited to study small systems, up to 4 x 4. 
For this reason, we also explain how one can go beyond that, by using the idea presented in [18]. The gauge conditions are, 
indeed, fixed in a RG inspired way, so to reduce the number of degrees of freedom. We can use a 4 — > 1 coarse graining scheme 
as follow. We cover the lattice with independent A s by chosing one of each two crosses and perform a RG step by keeping only 
the corresponding 6 eigenvectors with eigenvalue 1 of the corresponding A s , from the available 16 states. This reduces the 2 3 2 
degree of freedom of the 4x4 lattice to only 8 6 . The 8 remaining crosses still act on four coarsegrained sites. Now, we can 
repeat the procedure by fixing two further crosses and select among the 6 4 states the corresponding 6 x 3 4 . In this way, we can 
perform exact diagonalization on the two effective spins. The unused gauge constraints are projected at this level. This procedure 
is sketched in Fig. 7. If we want to deal with systems larger than 4 x 4 we cannot use the above construction. However, the 
introduction of an addtional RG step, which projects onto the low-energy subspace of the full Hamiltoniana, allows to study 
systems up to 10 x 10 (see, e.g., [53]). This can be done numerically as described in detail in [18]. 



III. QUANTUM SIMULATION OF THE U ( 1 ) GAUGE MAGNET 

In this section, we discuss the quantum simulation of Abelian, in particular U(l), gauge magnets in optical lattices. As the 
gauge magnets may be regarded as a prototype of lattice gauge theories, many of the major issues of the implementation are 
generic to any gauge theory. 



A. Many-body interaction through the Rydberg gate 

An important feature common to gauge theories on a generic lattice is the occurrence of interactions between several links 
together, involving at least as many links as the lattice coordination number. As such links appear in the dual lattice as sites, it 
implies the existence of a long range interaction between several sites at once. The aforementioned feature is very hard to be 
simulated. For instance, it is very hard to engineer in a controlled way such an interaction between cold atoms loaded in optical 
lattices. Such problem disfavors the analogic simulation of lattice gauge theory, i.e., the implementation of the gauge theory 
Hamiltonian as the effective Hamiltonian of the controlled system. 

The alternative is to consider a digital quantum simulation, i.e., the simulation of time-evolution of gauge theories via quantum 
gates. Its realization requires the capability of assigning a certain phase to a set of links,phases determined by the expectation 
value of each multilink operator entering the Hamiltonian of interest. The gates that implement these operations are generaliza- 
tions of C-NOT gates, and the key role of control over links at distance is played by Rydberg atoms [54]. 

Rydberg atoms are highly excited states that are long lived and have a huge dipolar moment [55, 56]. The latter gives raise 
to the phenomenon known as blockade: once a Rydberg atom has been laser induced, a volume of several jim is affected by 
its dipole field causing level shifts in the surrounding atoms, and, hence, forbidding the creation of another Rydberg excitation 
in such volume. Several experiments have demonstrated that it is possible to create entanglement over an extended region in a 
(laser) controlled way by exploiting such strong dipolar interaction [57]. 



14 




Figure 8: Schematic of the Rydberg gate. In the block of the left, the geometric disposition of the ensemble and control atoms is shown. The 
ensemble atoms, which admit as logical states |A) and \B), are placed sufficiently close to the control atom, whom logic states are |0) and 1 1), so 
that the Rydberg interaction, i.e., the dipole interaction in the Rydberg state of the control, is strong (shadow green volume). At the same time, 
the atoms are sufficiently separated to allow individual single-site addressing by laser pulses. In the block of the right, the pulses, needed to the 
conditional transfer from the logical state |A) to the logical state \B) of the ensemble atoms, and the level scheme are displayed, a) Sequence of 
the pulses (not in scale). The 7T-pulse of Rabi frequency Q. r (in red) is applied on the control atom only and induces resonant transition between 
level 1 1) and the Rydberg excited state \r). Such pulse is short and applied twice, at the beginning and at the end of conditional transfer. The 
pulse of Rabi frequency fl c (in blue) and the ^-pulse of Rabi frequency £l p (in cyan) act both on all the ensemble atoms simultaneously. The 
former is shining continuously and much more intense than the latter, which is smoothly ramped up and down, b) Level scheme for the control 
initialized in |0). In this case, the pulse fl r does not play any role. The 7T-pulse £l p connects with an (off-resonant) two-photon transition 
between the states |A) and \B) of the ensemble atoms. The frequency of the pulse Cl c is such that Q. c + Sl p connect the logic states |A) and |B) 
to the Rydberg excited state \R) with a resonant two phon transition. As Q. c 2> £l p , the Electromagnetic induced transparency scheme (EIT) is 
realized and the transfer between |A) and \B) inhibited, c) Level scheme for the control initialized in 1 1). In this case, the pulse £l r excites the 
control to the Rydberg state \r). Due to dipole interaction the Rydberg state of the ensemble \R) gets shifted in energy by V and it is no longer 
in a two-photon resonant transition with the logical state |A) and \B). Hence, the transfer between |A) and \B) occurs. 



Indeed, due to the blockade, the presence/absence of a Rydberg excitation plays the role of an ideal control qubit, as it may 
trigger or forbid laser stimulated transitions in all the atoms in the blockade volume, for instance all the atoms trapped in a unit 
cell of a lattice. The explicit realization of this idea has been pushed forward in [32], where the authors propose the mesoscopic 
Rydberg gate scheme, which is sketched in Fig. 8. Let us review briefly the main feature of this gate: for more details and for the 
analysis of the fidelity we refer the reader to the original reference. Let us consider a set of 1 qubits, one of which playing 
the role of the control. The logic states of the control {|0) , |1)}, and of the remaining ensemble qubits, {|A)|(., k = l,...,N, 
are hyperfine states of atoms. The goal is to achieve a protocol such that the ensemble states are untouched when the control is 
in |0) while induces a flip |A),t f-> IB)^ when the control is in 1 1). In formula, this reads 

G=|0)<0|®% + |1)<1|® rSKJ' (42) 

where Ijv is the identity in the space of the ensemble qubits and a x is the Pauli matrix acting as O^A)^ — \B)k, of \B)k = |A),t, 
\/k = 1,... ,N. We may suppose that both the control and the ensemble atoms are trapped in different sites of an optical lattice 
and that we have single site addressing. All the atoms can be excited by laser pulses to Rydberg states, which we indicate as \r) 
and \R) k for the control and the ensemble atoms, respectively. \r) and 1^)^. can be in principle different. The key idea is to use 
a Electromagnetic Induced Transparency (EIT) [58] in a modified Lambda scheme. Three different 7T-pulses are employed with 
Rabi frequency £l r , Q. c and Q. p , respectively. The first one is applied to the control at the beginning and at the end of the process 
and it is resonant between |1) and \r). This implies that the control is excited to the Rydberg state only when it is in the logic 
state 1 1). The other two lasers are such that the sum of their frequencies is exactly resonant with the transition \A) k , \B) k — > \R) k 
when the control atom is not in the Rydberg state |r), while it is off resonant when the strong dipole-dipole interaction is present 
(a more accurate discussion should take into account also the effect of simultaneous Rydberg excitations in the ensemble atom as 
done originally in [32]; such extra-dipole interactions can be kept under control and the functionality of the gate is not spoiled). 
Together with the condition A S> £l c ^> Cl p , where A is the smallest detuning of the p-laser from an intermediate atomic level 
\P), the above construction ensures that no transfer between \A) k and \B) k occurs when they are in two-photons resonance with 
\R) k , i.e., the control is in |0) ; on the contrary the transfer occurs with efficiency one in the absence of the two-photons resonance 
i.e., the control is in |1). This ideal result is not spoiled once the losses due to spontaneous decay of the long-live atomic states 
are taken into account. 

Before moving to the applications of such gate to the simulation of the gauge magnets, two comments are in order. First, the 
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Ensemble atoms 
Control atoms 




Figure 9: Schematic of the simulator. As indicated in the legend, the red dots represent the ensemble atoms, while the blue dots represent 
the control ones. Between the control atoms, we distinguish the ones at the lattice sites (intersections of continuous lines) and the ones at the 
center of the plaquettes. The former are needed to impose the gauge condition and interact with the ensemble atoms at the corners of the green 
squares. The latter control the dynamics of the plaquettes, which involves the ensemble atoms at the corner of the yellows squares. 



mesoscopic Rydberg gate can be modified to get any G' of the form 

G' = |0)<0|®% + |1)<1|® ((ZsA*^* ) =6dU k (fi k ) (|0)(0|®% + |1)(1|® I ®* r ?°i' I )(g)C/J(%), (43) 





where a s is any real number between and 1%, n k a generic 3-vector of modulus 1, G = (c7 x ,a y ,<J z ), and Uk{nk) the unitary 
transformations in each single qubit space. Indeed, the a ? dependence is simple obtained by considering an a ? -pulse for the 
lasers c and p instead of the 7£-pulse considered originally in the mesoscopic Rydberg gate. On the other hand, in order to get 
the a combination n k ■ G, the single site addressability is exploited to perform the appropriate rotation of each qubit basis. This 
operation obviously amounts to apply single qubit gates on the ensemble atoms before and after the G gate. 

Second, as originally shown in [33], the mesoscopic Rydberg gate can be used to read out the expectation of a (two-valued) 
operator on the ensemble atoms and transfer the information to the control qubit, once combined with rotations on the control 
itself. Indeed, let us consider the effect of the unitary transformation 

G = exp[/|^]Gexp[-/|^'], (44) 



once the control qubit is initialized in |0). For a generic state of the ensemble atoms |A) = Of s |A+) + y 1 — a^|A_), where 

with \X±) are eigenvectors of O e = ®£<7?, O e \X±) = ±\X±), and \a g \ < 1, the final state of the control is computed to be 

|c) = <X g \0) + y/l — ce||l). Furthermore, by inverting G, any single qubit operation applied to the control can be mapped to 

the ensemble state, i.e., on the logical qubit composed by |A + ) and In case we are interested to read a more complicated 
operator, as in (43), appropriate rotations of the ensemble qubits are required as well, & — exp[/|a^] G' exp[— z'f Cc]- 



B. Dissipative preparation as ground-state preparation of frustration-free Hamiltonian 

We are now ready to discuss in full detail the protocol for the simulation of the gauge magnets. The link qubits are identified 
with the logical states of the ensemble atoms (see Fig. 9). For convenience, we choose the following relation between the 
{ | +) , | — ) } link basis of section II C and the logical (atomic level) basis { \A) , \B) } 

l ± ) = ^>. (45) 

The first conceptual and practical step is to impose the gauge condition (7), i.e., to reduce the possible states of the system to the 
gauge invariant ones. It is worth to notice that this is equivalent to find the (highly degenerate) ground-state of the Hamiltonian 
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Hcond = — \Hs{As(g) (,?))> with <? a generic element of Z^v or of U(l) such that g" ^ e for n <N,oi for any n, respectively. 
In the U(l) case, we may choose H conc ( to be 

1 . CCg , CCg y . CCg y ■ \ 

H cond = £ /i^ om/ = - - £ e' T CT 'i ® e' T CT '2 ® e - f °? 3 ® e - x °? 4 +H.c, (46) 

s Z s 

where a ? is any integer, for instance 1. 

Independently of the specific election for H con d, it is, evidently, a frustration-free Hamiltonian: each local term commutes 
with the others, and the ground-state can be obtained by minimizing each local term separately. In this situation, a dissipative 
preparation of the ground-state, as originally proposed in [38, 40], and applied to quantum simulation in [33] (for a very recent 
review see [59]), is appropriate. The idea of dissipative preparation is very simple and corresponds to a classical cooling. That 
is to say that each time the local state has no minimal energy is transformed to a minimal energy eigenstate by a local operation. 
Let us detail the procedure in the case of (46). Our aim is to read out the expectation of A s {h s cond ) on the local state and to 
change it, for instance by flipping an ensemble qubit, if is not satisfying the gauge condition (or equivalently (ft L_j) = — 1). As 
our local state is given by the four links attached to a site and A s is a unitary operation that can be conditionally implemented 
(see discussion above) by the mesoscopic Rydberg gate, the required operation is obtained with the gate sequence 

U D = G(A s fU^G(A s ), (47) 

where G(A S ) = e'f °? (|0)(0| ® 1, + 1 1)(1| ® A^e' 1 * = ±l c ® (l s +A S ) + ±CT* ® (A s - t s ), and 

Uf' 9 = \0)(0\®± s + \l)(l\®Z((p), 

is the conditional flipping of one of the links entering s with probability <p, let say s\, £(<p) = e ,(pfIs i . 

Two comments are in order. First, we note that the spin flip E(<p) is not necessary projecting the local state \y/) s to the 
subspace (A s ) = 1 in one step, even for <p = n/2. Furthermore, if h s cond \y) s = — cos(2a ? )|y/). s , it follows that h s C(md L(^ )\y/) s = 
— cos(a g )E(|)|y/) 4 but if h s cond \y) s — — cos(a ? )|y/) 4 -, it is found T.(j)\y/) s has component on the —1 and — cos(2a g ) subspace 

with amplitudes y| and j [60]. That is to say, the "energy" is not monotonically decreasing at each step. However, the relation 
above is telling us that, by iterating the process, the probability of the state to be excited is exponentially suppressed, as the 
probability of being reheated scales like {\) n l 2 , n = # of steps. 

Second, it is worth to notice that the operation Ud is unitary, but the dissipation enters in the game as, in order to iterate Ud, 
we need to reset the control qubit to |0) c . As explained in [33], such decoupling of the control dynamics can be obtained by 
inducing a radiative decay of |1) ( .. This implies that the final result of one dissipative iteration is a mixed state. Indeed, by 
naming \ y/) s any initial pure local state and assuming that the control qubit is initialized in |0), the action of Up is 

|0> ® W) s ^ |0> ® |0o).v + 11> ® |fr>„ (48) 



where 



This means that, in consequence of the decay of 1 1) to |0), the density matrix after a dissipation step evolves as 

\0)®\V)(Y\ s DiSSiPati(m > IftXftl + lftXfrlr (3D 

The most efficient cooling procedure is for (p — 5. However, as discussed in [33], it is possible to show that, in the limit of small 
<p, the above cooling procedure is governed by a master equation of the Lindblad form. By defining the jump operators 



(t-At )af(l-A s ) t (1+4K(1-A S ) 
"4 =c 0' c i = 4- 



c = =c , ci = , (52) 



the expansion of (49) and (50) in <p around <p = takes the form 

,2 



\(j)o) s = ^l-/(pco-y(cg + c| Cl ) + <9((j!> 3 )^ \y) s , \<j)i) s = (-i<pci+0(<p 2 ))\\i/) Sl 



(53) 
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which implies that the variation of p after one cooling step can be written up to second order in <p as 

Ap=Atp = -i(p[c ,p] + (p 2 (c Q pc Q + Cl pc\ - ^{(co+c{ Cl ),p}j +0(<p 3 ), (54) 

where At is the time needed for a the single time step. It is immediate to see that the above equation has no other fixed points 
than a density matrix built up of pure states, p = |V0')(V0'l> satisfying A s \ Yj) = \ In comparison to [33], we find that 

two jump operators are associated to the cooling procedure and the presence of a coherent term linear in <p, which is not affecting 
the convergence as it is zero on the gauge invariant subspace. 



C. Unitary evolution of the U( 1 ) gauge magnet 

Once the gauge condition has been successfully implemented, we can focus on the dynamics of U(l) gauge magnet, governed 
by the Hamiltonian (18). As explained above, our goal is to implement the unitary evolution U t = e~' Ht in terms of local unitary 
gates [33]. The simplest way to do it is to exploit the Trotter expansion to write 

Ut-Ylie-^e-^"^)* , (55) 
v 

where H p and Hj are the plaquette, H p = — cos 9(B p +B y p ), and the link//; = sin 6(X/ +xj) terms, respectively. As the Rydberg 
gate are sufficiently fast (Af ~ fis) compared with temperatures achieved for the atoms in the lattice, the above approximation is 
consistent and we can treat each term as commuting in the time-evolution. 

Let us focus first in the implementation of the link part. As the X[ commute between them, such term is just a product of the 
single qubit rotation U\ = e -'(i+eosa g )r e -i(i-cosa^)rc7/ applied to each link. Thus, its implementation is straightforward. 

On the other hand, the plaquette term is slightly more involved. The easiest way to treat this term is to expand explicitly the 
exponential into a sum of tensor product of a. In order to reduce the number of terms to be simulated, we may use the Trotter 
expansion further. As Z and can be written as Z = a ^ 2 ' a an d Z T = a ~t' CT , respectively, we have 



cos 9 £ U) 



h p = — 3— LQp 

° 7=1 

-((J ^ )® 2 0(a- v )® 2 -((J y )® 2 (8)(a z ) 8,2 ), (56) 

which gives 

e-'^^fle'^Qp^. (57) 

7=1 

Now, each of the above unitaries can be implemented using Rydberg gates as explained in [33], by applying the sequence 

js&dP* = G(Q^e'^^G(Q { p j) ), (58) 

for G(Q p j) ) = expfif cr e y ](|0)(0| ®t p + |1>(1| ®Q { J ) )ex^[-i^ a y c ], which is of the form of (43). 

It is worth to conclude with a final remark on the preservation of gauge invariance during the above process. In principle, as 
the Hamiltonian commutes with the gauge condition (7), the time-evolution conserves the gauge condition, i.e., it is not driving 
the system out of the gauge invariance subspace. In fact, in order to minimizing possible violations of the gauge condition due 
to dephasing it is worth to keep applying the dissipative process while other unitary operations are performed. By doing so, the 
unitary evolution is not altered. 



D. Proposal for the preparation of the U( 1 ) gauge magnet ground-state 



We are finally ready to discuss the preparation of the U(l) gauge magnet ground-state. As in the link dominated phase the 
ground-state is a fully polarized product state that can be efficiently prepared using single site addressing, we focus on the 
plaquette dominated phase. Our plan is to realize in practice the procedure explained in section II C. Hence, the idea is to 



18 



construct the ground-state in two steps. First, we apply the dissipative procedure to the frustration-free Hamiltonian obtained 
considering half of the plaquette (for instance the yellow ones in Fig. 4) in order to find the corresponding ground-state |£2(0)) 
defined by (29) and (30), and, then, we drive such state, adiabatically or using the CRAB procedure [42], to the ground-state of 
the full plaquette Hamiltonian. It is worth to notice that prior to the ground-state preparation it is not necessary to impose the 
gauge condition as |£2(0)) is gauge invariant by construction. However, for the driving, the final remark of the previous section 
applies. 

The procedure explained above works for any chosen background of static charged matter, once the proper initial state is 
prepared and the right gauge condition is applied at site. As in section II C, we focus on 1) the preparation of the ground-state in 
the absence of charges and on 2) the preparation of the ground-state in presence of two isolated + and - charges. 

1. Preparation of the ground- state in the absence of charges 

The preparation of the initial state to be driven can be obtained by applying the discussion of section IIIB to the local 
Hamiltonian H p , p ejyellow plaquettes}. We note that the unitary operator Uh p = exp[z'§ (B p +Bt)] acts on the eigenvector of 
H p (see definitions (30-32)) as 

U Hp \l) P = \l) P , U Hp \-l) p = -\-l) p , U Hp \0) p = -i\0) p , (59) 

i.e., in a very similar way to the star operator A s . That is to say, we can get the right dissipative process just by substituting A s 
with Uh p in the gate derived in section III B. However, the C-NOT gate |0)(0| ® t p + | 1)(1 1 <8> Uh p is hard to be implemented 
exactly (i.e., without applying Trotter approximation) using Rydberg gates. 

We can bypass easily such difficulty by exploiting the exact knowledge of |£2(0)) = <8>pG{Yeiipi } \^) P < where the local state 
| \) p can be explicitly written by ordering the links of plaquette anticlockwise starting from the bottom as 

|l) P = -^(|+,+,- ,-), + |- -,+,+),)• (60) 

We observe that such a state is completely determined by: i) being a combination of states with an equal number of + and — , ii) 
being the stabilizer state of both O x = O x ® O x ® 1 ® 1 and O z = O z ® O z ® O z <g> O z i.e., satisfying 0^\\) p — \ \) p , fi =x,z. 

The condition i) can be realized in different ways. As suggested in section II C, the most direct one is to prepare the product 
state |y/ ) = G 

{Yeil.pl.} I~h ~K — ) — ) P hy exploiting single addressing. As an alternative, we could first impose the gauge 
condition and, then, repeat the dissipative procedure of section III B by choosing as local state the (yellow) plaquettes, as unitary 
operator A p = expf/^a*]® 4 , and as flip operator the simultaneous rotation of two plaquette's qubit by a z , £ = e'f ' ® e I¥>0 2. 
Indeed, the gauge condition ensures that the initial plaquette state has an even number of + and the double flip converts states 
with zero o four + in a desired state. 

The condition ii), once i) is realized, can be achieved by applying the Kraus maps for Bell state pumping engineered in [61]. 
Indeed, let us suppose that we have prepared the state |y/b)- Hence, we have to impose by dissipation only (O z ) = 1 on the 
plaquette state as O x \y/o) = 1 by construction. The dissipative gate 

G(0Z)ttfz,9=§G(0Z) + radiative decay] 1) ->■ |0) 

does the job in one step. 

Now, we are ready to turn on the remaining plaquette interactions. This can be done by increasing the coupling constant c w 
in (41) adiabatically from to 1. From the digital simulation point of view, this means that unitary time-evolution procedure of 
section III C has applied to the plaquette Hamiltonian (41), where c w is changed smoothly, for instance linearly, from an iteration 
to an other. However, as first shown in [50] and reviewed in section II C, in the infinite dimensional system, the plaquette 
dominated phase becomes gapless as c w approaches 1, with the gap that scales as \. This means that the adiabatic driving for 
a finite system is still possible but it becomes harder and harder as the dimension of the lattice increases. This implies that the 
ramping time has to be taken bigger and bigger. In order to reduce such time, optimal control techniques can be used. Between 
them, the so call CRAB method [41], which has been successfully employed, for instance, to dramatically reduce the time 
needed to bring a superfluid gas into a Mott insulator state [42], seems very promising as it allows for closed-loop optimization 
experiments. The idea is to find the optimal function c w (t), satisfying c,„(0) = 0, such that after a time T minimizes a cost 
function which in our case is a linear combination of the energy of the driven state computed with the final Hamiltonian 
c w = 1 and the fluences, a measure of the oscillation of c w {t). The novel feature of the CRAB method is that the optimal driving 

at the step /, c„ (?), is obtained by minimizing with respect to few free parameters in the driving probes c[' v (t) £ Pkgkif). 
The gk(t) are element of a favorite orthonormal basis, for instance the Fourier basis gk{t) = cos(fc|if ). In order to further improve 
the convergence, the frequencies Wk = kj may be randomly perturbed Wk — > Wk + i~k- One then feeds the result of the optimization 
back to the algorithm so that the final energy is improved further up to convergence. 
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Figure 10: Comparison between linear and CRAB protocol drivings. The preparation of the ground-state of the gauge magnet Hamiltonian 
(41) is performed in a given time, shorter than the adiabatic one. The blue line represents a linear ramp from c w = to c w = 1. The red curve 
represents the result of one optimization cycle of a CRAB preparation with seven random frequencies as free parameters [41]. The horizontal 
lines represent the energies of the ground-state and the first excited state. It can be appreciated how the linear map fails completely to deliver 
the desired ground-state while the CRAB map gets very closed to it. 



We have numerically tested the CRAB method up to a 4 x 4 lattice, where we are still able to compute the exact ground- 
state for c w = 1 using exact diagonalization. One has to keep in mind that the energy of the final state of the time-evolution 
is monotonically decreasing with the number of iterations. Since each iteration takes very long (typically around a week on a 
standard PC) we have decided to stop our optimization algorithm as soon as the final energy got well below the energy of the 
first excited state. This typically happens after only one optimization cycle. For the sake of comparison we report in Fig. 10 the 
energy of the state along its time-evolution and compare it with respect to the energy of the state obtained using a linear ramp 
that goes from c w = to c w = 1, in the same amount of time. It can be clearly seen how out of the adiabatic regime the linear 
ramp fails completely to prepare the ground-state, as the evolution is by far too fast, while the CRAB optimization with only one 
iteration gets very closed to it. 



2. Preparation of the ground-state in presence of two isolated + and - charges 

In presence of two isolated + and - charges, two main difference in the preparation occur. First, the initial state to be driven, 
i.e. the ground-state of the system with only half of the plaquettes turned on (the yellow sublattice of Fig. 4), is different. 
Indeed, as explained in section II C, it is again a product state of eigenvectors of the plaquettes turned on, but it differs from the 
ground-state in the absence of charges on few plaquettes, which we c ailed flipped. On the flipped plaquettes, gauge invariance 
forces to take an eigenvector of eigenvalue 0. Furthermore, the flipped plaquettes can only occupy a continuous segment of 
a southeast-northwest diagonal of the sublattice, see Fig. 5. Accordingly, the charges can be placed consistently with gauge 
invariance only on the sites of the first and last plaquette of the segment. To be precise, a charge can occupy all the sites except 
the one shared with the contiguous flipped plaquette of the interior. Hence, by definition the flipped plaquettes of the interior 
contains no charges. The position of the charges fully determines which are states on the flipped plaquettes. Let assume for 
convenience that the positive charge is on the first flipped plaquette at the bottom-left and the negative one is on the last at the 
top-right. It follows that the proper state for the plaquettes in the interior of the segment is | + + ++). For the plaquettes on the 
frontier we have two possibility: the charge occupies the site, d, on the diagonal of the segment (not shared with other flipped 
plaquettes), or it occupies one of the two off diagonal. In the former case the state is again | + + ++). In the latter, the link 
connecting d with the site of the charge is in |— ) with all the other links are in |+). Hence, it follows that the plaquette at the 
bottom-left is in the state | — if the charge is below the diagonal or in |+ + +— ) if the charge is above. Similarly, the 

plaquette at the top-right is in the state |H h+) if the charge is below the diagonal or in |+H h) if the charge is above. 

Such states can be easily prepared by single-site addressing or using dissipation. It is worth to notice that the configuration with 
the positive charge at the top-right can be obtained from the one above by combining a ^-rotation of the lattice, which puts the 
charges in the desired configuration but reverses the orientation, with the flip of the state of each link, |+) f-> |— ), which restores 
the right orientation. Hence, the preparation of initial state to be driven in presence of charges does not imply much more effort 
than the preparation of the not-charged one. 

The other difference in preparation of the ground-state with charge with respect to case without is that the gauge condition 
to be enforced during the driving in order to minimize the error is different. Indeed, on the site where ±1 charge is placed the 
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gauge condition reads 

A s± \ W )=e ±,a \ W ). (61) 
This means that the dissipative procedure has to force the eigenvalue 1 of the operator A' — e f,a A s± . By noticing that 

e ±''f ^ G(A) = e^GiA') = e ±; f (|0)(0| ® 1 + |1)(1 1 ® e T '°A) , 
we conclude that the correct gauge condition can be imposed via the modified version of dissipative sequence (47) 

U±=G(A>JU^G(A>J, (62) 

where G(A ! S± ) = e't^ e ±! 'f CT c : G(A 1 . ± )e"'i°^ = e ±! '" ML e ® 1,;± +e J Aj± + og eT ^ , as the phase e ±! '? cancels out in 

the gate. 

E. Validation and usefulness of the simulator: absence of symmetry breaking and confinement 

Once the proposed experiment has been realized, one has to be able to ensure that the simulator is indeed behaving as expected, 
and thus validate the results. This is in general a non-trivial task [62]. In the case of a gauge theory simulation, there are, however, 
at least two expected properties that one could use as a "smoking gun" of the desired physics. 

First, one can check that the expectation value of any non-gauge invariant operator is zero as predicted by gauge invariance for 
well defined, i.e., non anomalous gauge theories like the gauge magnet under scrutiny. The simplest of such operators are (twice) 
the Hermitian and antiHermitian part of Z, in the logical basis, Z + Z^ = a z and —i{Z—Z^) = —a y , respectively. Verifying that 
the expectation value of these operators is zero for any link provides a simple necessary condition of the consistent realization 
of the gauge symmetry - obviously, it is saying nothing about how good is the simulation of the dynamics. In order to make 
the check of gauge-invariance more compelling - it is well known that while the expectation value of local operators can be 
reasonably close to what craved, the state can still be orthogonal to the one desired- we may consider similar operators on 

larger support such as non-intersecting string of links, ®/{T, . The vanishing of all the (0;a, ) would provide a stronger 
confirmation about the gauge invariance of the prepared state. 

It is worth to notice that, as these measurements are destructive, i.e., they are non-commuting with the gauge-condition and 
with the Hamiltonian, the (ground-)state of the model has to prepare again after any of these measurements. 

The second important check is non destructive, and concerns the dynamics of gauge theories. As argued in section II C, a 
common mark to them is charge confinement. This can be tested by measuring the charge gap defined as the difference in energy 
between the ground-state with opposite ±1 charges and the vacuum. It requires the comparison of the ground-state energies 
in presence of ±1 static background charges at different distances. The vacuum may be thought as the limit configuration of 
zero displacement between the charges. As explained in section II C, the different ground-state preparation amount in changing 
the gauge symmetry condition from (7) to (34) at the two sites where the charges are located. In presence of confinement, the 
dependence of the charge gap upon the inter-charge distance (for the model proposed) has to be linear, with the proportionality 
factor expressing the string tension, see (36). This measurement can be performed at both stages of the ground-state preparation. 
Indeed, after the dissipative stage, charges should already be confined and the string tension in known to be a = 1 /V2. The 
value extracted from the experiment can, then, be used to check whether the first step of the state preparation has been successful. 
After the implementation of the ground-state of the full Hamiltonian, for instance, via the CRAB optimization, charges should 
still be confined but with a different and unknown value of the string tension. At present, this value cannot been determined with 
other means for sufficiently large lattices (the string tension itself is not well-defined for small lattices). Hence, it would be first 
interesting determination of the quantum simulation. 

F. Some remarks on the experimental set-up and summary of the experimental sequence 

For sake of clarity, we conclude this section by summarizing which are the key ingredients to build up the quantum simulator 
and which is the experimental schedule of the simulation. 

1. Requirements of the implementation 

• The first requirement for setting up the experiments is the choice of ultracold atoms suitable for the implementation of 
the mesoscopic Rydberg gate and for the load in the appropriate optical lattice. In the simplest scenario we may consider, 
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ensemble and control atoms are of the same atomic species, for instances Rubidium. As a consequence, an ordinary square 
lattice, obtained by shining the atoms with two pairs of counterpropagating lasers in the X and Y direction, respectively, 
is sufficient to host all the atoms. If distinct atomic species for ensemble and control are considered, a more complicated 
superlattice structure has to be employed, as, for separate, physical (link qubits represented by the ensemble atoms) and 
ancillary (control atoms) degrees of freedom occupy two square lattices at 45° (see Fig. 9). In any case, the lattice 
potential has to be sufficiently high that the free hopping is suppressed (compared with the time scale of the experiment): 
for simplicity, the atoms are assumed to be in Mott state with just one atom per site. 

• The second crucial requirement is single-site addressing, i.e., the capacity of (laser) manipulating the atom in each site 
individually. Single-site addressing is technically hard but possible [63] in ultracold atom experiments. It is at the hearth 
of the functioning of mesoscopic Rydberg gate, see Fig. 8. Furthermore, due to single-site addressing, the position in the 
lattice is sufficient by itself to distinguish control atoms from the ensemble ones. 

• The third requirement is that the system be sufficiently cold such that the energy scale E associated to the Hamiltonian 
can be resolved, i.e., KT < E where T is the temperature and K the Boltzmann constant. The energy scale E, roughly 
speaking the normalization of the plaquette term (for convenient it has been fixed to 1 in our analysis), is limited by 
the Trotter approximation to be much less than the inverse of the time At in which one unitary step of time evolution is 
performed. Such time At is determined by the number of sequential Rydberg gates needed to engineer the plaquette 
term evolution, At = NrIr , as the further delay due to single-qubit rotations entering the process is negligible compared to 
the Rydberg gate delay f«. Hence, it follows A^ <C tgKT. For state of the art experiments t^KT > 10 2 , as tR is about few jxs 
and T is about few tens of nK. Let us now compute Nr. The plaquette term decomposes in eight elementary monomials of 
a 11 , jX = 0, . . . , 3, (J° = 1, as shown in (56). Two of them, (c z )® 4 and (a >- )® 4 , are commuting when applied on contiguous 
plaquettes. The other six, (a z <g> ff- v )® 2 , (a y (g) a z )® 2 , o z <g> o y <g> a y <g> a z , a y ® a z <g> a z <g> a y , (ct z )® 2 <g> (a- v )® 2 , (a y <g> o z )® 2 
can be decomposed in twelve commuting monomials once applied on the alternated plaquettes. As the unitary evolution 
of each monomial involves two Rydberg gates, the total number of gates performed is Ng = (2 + 6 X 2) X 2 = 28, which 
is sufficiently small to satisfy the condition. Performing gate optimization, A^ may be further reduced. 

2. Experimental sequence 

1. Loading of the ultracold atoms in the appropriate square (super)lattice as in Fig. 9. 

2. Preparation of the initial appropriate product state as in section III D, by using dissipation via Rydberg gates (or dissipation 
combined with single-site manipulation of individual ensemble atoms). As discussed in section II C, the matter content of 
the final desired state to achieve (after the driving, see next point) determines which is the appropriate initial state to start 
with. 

3. Adiabatic unitary evolution using Rydberg gates as explained in section HID. Commuting operations can be performed 
on parallel. More efficient driving as the CRAB method can be adopted. 

4. Error minimization by enforcing the gauge condition during the driving via dissipation (see section III B). The two proce- 
dures, driving and dissipation, can be applied for alternative intervals of time. 

5. Validation of the simulator and "computation" of the simulation. Once the ground-state of the system (for the chosen 
matter content) has been achieved, it is characterized by measuring correlation functions. Such operation can be performed 
by mapping the expectation value of the operator on the control atoms with the gate G, as defined in (44). Indeed, 
the relevant operators to be evaluated are (linear combination of) monomials of <7^, Q = G^ 1 ® ■ ■ ■ <X> (J^ k , both for the 
validation and the evaluation of the string tension. The latter requires the evaluation of the energy decomposed in the 
sum of commuting monomials with k = 4 (as explained above). In practice, one starts with the control initialized in 
|0). Once applied the unitary map G, the control will be in the linear superposition |0) and |1), reflecting which is the 
decomposition of the state in eigenvalues ±1 of Q (see explanation below (44)). The control state can be measured in 
various ways, for instance using florescence [64]. Afterwards, the original state of the simulator has to be prepared again, 
as the measurement alters it. In order to get enough statistics, the whole process has to be repeated several times. 

IV. CONCLUSIONS AND OUTLOOK 

In this work, we have introduced a new derivation of the U(l) gauge magnet in 2+1 dimension as a special case of a generic 
Abelian LTG, and discussed its phase diagram for a natural choice of the gauge invariant Hamiltonian. The advantage of such 
derivation is that it shows explicitly that the U(l) gauge magnet may be viewed as a truncated formulation of the standard 
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Kogut-Susskind U(l) LGT theory. In particular, the amount of resources needed for the simulation is drastically reduced as the 
local Hilbert space shrinks from infinity to 2D (qubit). This allows to engineer a digital simulation of gauge magnet in optical 
lattices based on the latest developments in Rydberg gates. 

In this context, we have proposed a protocol that allows for both preparing the ground-state of the model and performing 
its out of equilibrium dynamics. We also discuss two possible measurements that can provide sine qua no conditions to the 
correctness of simulation. In order to perform a successful simulation, one has at least check that the states obtained are gauge 
invariant and that static charges are confined. The determination of the string tension for the full-fledge Hamiltonian of the 1/(1) 
gauge magnet would also provide the first non-trivial result obtained from a quantum simulator in the area of gauge theories. 

Our results pose interesting questions and open novel possibilities that deserve further studies. First of all, we (re)discover that 
the simplest plaquette Hamiltonian gives rise to confined phase, as expected, but gapless. Such property distinguishes the gauge 
magnets from the ordinary U(l) gauge theory in 2+1 (cf. [52]), and it is a consequence of novel gauge condition arising from the 
different choice of the Hilbert space (such choice, for instance, rules out the existence of short closed-string configurations). How 
the choice of the group representation affects universal properties like phase diagrams and phase transitions is a fundamental 
question that raises naturally here and is asking for solutions. 

Similarly, our constructive approach to gauge magnets allows an infinite class of Hamiltonians built up of local terms to be 
considered. It is totally unknown, at the moment, whether novel phases and phase transitions are described by such models, and, 
for instance, whether the Polyakov phase could coexist in the same phase diagram with the gapless phase we have encountered 
in our study. The above questions are especially stimulating as nowadays novel tools in both classical simulation algorithms and 
in quantum simulation may be designed and employed to find a solution. 

Furthermore, another interesting line of investigation is the search for analogues of the gauge magnets formulations in 3+1 
dimensions and/or for non-Abelian groups. In parallel, the introduction of the coupling to charged matter is very appealing. As 
relativistic fermions are relatively easy to be simulated analogically in optical lattices, a combined digital/analogic simulation 
approach is desirable. 
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Appendix A: Gauge group Z2. 
1. The local Hilbert space 

The Z2 is the simplest LGT one could imagine (apart from the percolationg LGT of [66]) where the orientation of the lattice 
is not important, since there are only two different elements of the group, the first is 1 and the other is e, which has the property 
that e 1 = ee^ 1 = 1. The local Hilbert space with spin 1/2 is the algebra of the group and the representation matrices of the group 
coincide with the regular representation of the group. The regular representation is defined to be the representation in which the 
matrices are constructed directly from the multiplication table of the group. In order to get the representation of multiplication 
by a given element, one substitute the given element in the multiplication with one and all the others with zeros. For the Z2 case 
we have the following group multiplication table 





1 e 


1 


1 e 


e = e 


e 1 



thus X(l) = 1, X(e) = o x . There is an important theorem relating the regular representation with the irreducible representations. 
This theorem states that the regular representation contains all the irreducible representations a number of time equal to their 
dimensions [48]. 

In this case, we have two A s (g), A s (l) = t and A s (e) = of 4 . Thus the condition (9) is projecting out half of the 2 4 states 
associated to a site. 
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2. The gauge invariant Hilbert space 

In order to show that the condition (9) is neither trivial or empty, and how it can be implemented, we first consider the minimal 
lattice, made of only one plaquette. This is shown explicitly in Fig. 11 v). The gauge invariant Hilbert space is embedded in the 
2 4 = 16 dimensional Hilbert space {1/1/2/3/4)}- Due to the periodic identification of of the links on the plaquette, (9) applied for 
si implies 

at®a} c %l 4 ) = \hk). (A2) 

Since a x has eigenvalues ±1. We have two possibilities I/1Z4) = |++), I/1/4) = | ), where 

05c|+) = |+), OEcl— > = — I— >- (A3) 

By iteratively applying all the A s (g) on subsequent sites, one realizes that there are two allowed states and thus a generic gauge 
invariant state is a linear combination of them 

|0+> = !+ + ++>, |0-) = I )• (A4) 

This formalize the naive intuition (which has to be modified for finite systems with periodic boundary conditions [18]) that gauge 
invariant states are related to the elementary plaquettes of the lattice. 



3. The gauge invariant operators 

The operators compatible with gauge invariance constraints are either 

• product of a x on arbitrary links of the lattice 

• product of Oj-t x on closed paths 

The simplest choice is then a[ l ® a' 2 ® C z ' 3 <g> o[ A . Going back to the lattice formed by a single plaquette of Fig. 1 1 v) we define 
the Hamiltonian 

H = -B P = -ai [ ^ a' 2 ® a, 1 ® <x! 4 , (A5) 

that is the form that the generic Hamiltonian of (18) on the chosen lattice when 0=0. In this simple example, we see that the 
ground-state is given by the linear combination 

IVo"> = ^ ! (l+ + ++} + l », (A6) 

and the first excitation is given by 

lvf P ) = -^(l+ + ++)H »• (A7) 

We also see that the two states are separated by an energy gap (that does not depend on the system size), 

A 2 = 2. (A8) 
On the other hand, if we consider the opposit limit of the Hamiltonian of (18) and set — 7C/2, we obtain 

Hn/2 = -L°i')> (A9) 

k 

which ground-state becomes 

IVo*> = 1 + + ++), (A10) 

and which first excited state is 

|VD = I ). (All) 

Already in this simple example we can appreciate that in the two limit considered of (18) we are in a gaped phase but the 
ground-states of the Hamiltonian in (A9) and (A5) have different entanglement properties. Indeed, the ground-state of the latter 
Hamiltonian is a product state, while the one of the former is locally maximally entangled, since the reduced density matrix of 
a single link is proportional to the identity. This reasoning can be extended to larger systems since both Hamiltonians of (A9) 
and (A5) are fixed point of the RG flow [18]. The general Hamiltonian (18) is not in general a fixed point of the RG flow. From 
the above discussion it is easy to accept the appearence of a phase transition at a given C in between 0=0 and = k/2 that 
separate the domain of attraction of the RG fixed point at = usually called the topological-deconfined phase with the one of 
the RG fixed point at = n/2 usually called the confined phase. 
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Appendix B: The case of the group Z3, prototype for the generic Z# 



i) 



L x L 



ii) 



Ml £ Vj 



t J f 



iii) 



•a 



A, (a) 



«(s)l<Wi, 



iv) 



v) 



SlJU siJX 1 \ I, I 



- ^ *-»-S2 

I ! 1 ^2 

'3 



_ li9 

A.(fl) 



Figure 11: i) Example of an oriented lattice Jz? with L sites .s and N = 2L links /. ii) A many-body state of the lattice is obtained by defining 
a local Hilbert space "V on each of the links. A state of the lattice is then a state of the tensor product \y) £ Mq = 'f® N . iii) For each site 
s of the lattice we can define operators A s (g) acting on the local Hilbert space 3V G \ S = y® 4 , where f is the Hilbert space of a link ending 
in i. A s (g) multiplies each link state \hi.) s for the X(g) or X(g~ l ) depending on the fact if the link is entering or exiting the particular site. 
Gauge invariance is the requirement that the state A s (g) |fy)s) = ®i l^/ ; )s))-i v ) A particular case is obtained when the local Hilbert space 
y is chosen to be the group algebra ^(G). In this case, the states are labeled by \gi,) s and changing the orientation of the link, corresponds to 
\hi t ) s — > \h i^s. v) A simple example of a one plaquette lattice, discussed in the main text, is worked out in details. It consists of four sites, 
si ■ ■ ■ S4 and four links / 1 • • • I4 . The four elementary gauge transformations related to the four sites are explicitly introduced. 



The next step of complication is obtained by considering a Z3. The main difference with respect to the Z2 theory is that the 
group algebra of Z3 has dimension 3 - N for a generic Z^ - and not all the elements of the group are the inverse of themselves. 
Hence, the orientation of the lattice plays a central role. From now on, we will consider two different basis, the group algebra 
basis 1 and the basis in which the regular representation becomes the direct sum of the irreducible representations \t). 

In the |g)-basis the X(g)s read 





1 


e2 = e x l 


e\ =e^_ 1 


1 


1 


ei 


e\ 


e\ 


e\ 


1 


e 2 


ei 


e2 




1 



(Bl) 
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(B2) 
(B3) 

(B4) 

As done previously for a x , we can diagonalize all the X and we pass to the basis |f,)(f, |, i = 1 . . .3. In this basis 

X(l)=diag{l,l,l}, (B5) 

X{e l )=diag{l,e¥ i ,e 2 r i }, (B6) 

X{e 2 )^diag{l,e^ i ,e'^ i }. (B7) 
We then can find the gauge invariant states in the |f) -basis 

\h.h,h,h), |?2,?2,?2,?2), |f3,f3,f3,*3)- (B8) 

Now, let us focus on the construction of the plaquette of the Hamiltonian, taking properly in account the orientation of the 
links. In order to do so, we have to specify the form of the Z operator for the group Z3, which generalizes the a z employed for 
Z2. In the |g)-basis, it is 

1 

cT j . (B9) 

,0 ef' 

It has the following properties 

ZX(e\)=e7 i X{e\)Z, ZX(e 2 ) = e^ i X(e 2 )Z. (BIO) 

If we orient the plaquettes in such a way that they always have one entering and one exiting link for each site (either anti- 
clockwise or clockwise) the plaquette will meet a gauge transformation on two consecutive links but the orientation is such that 
the gauge transformation on these links will look like X(e) ®X(e~ l ) . Since the orientation of the plaquette is different from 
the orientation of the gauge transformation given in Fig. 1 1 iii) we are multiplying one of the two links for X(e~ l ) instead than 
X(e). Hence, the only possibility to obtain an operator commuting with all the gauge transformations is to take again a product 
of Z as 

[Z^Z,X(e)^X(e- 1 )} = 0, (Bll) 

accordingly to (BIO). It follows that the correspondent Hamiltonian term, written on the anti -clockwise or clockwise oriented 
plaquettes (as the one of Fig. 1 1 v) ), reads 

H r = AL{ Z h z k z h z k + z i z ;A z i) • (B12) 
z p 

This is, as expected, the explicit realization of Hamiltonian of (18) with 9 = for the group Z3. It is important to notice that 
now the l\ ■ ■ ■ In are the links around a given plaquette in the order of Fig. 1 1 v). Once we move back to the standard orientation 
of the 2D space of Fig. 1 1 ii) we recover the standard form [47] 

R p = -lL( z i z i z h zi * + z 'i z '2 z i z i) ■ < m3 > 

Once the basis is rotated from \g) to |f), the form of the Z operator is determined by (BIO) to be 

Z:\ti)->\ti +l ), (B14) 

and 

tf:\t t )->\tt-l), (B15) 
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i.e.,X(ei) ->Z, and Z->X(ei). 

Again we can consider the ground-state of the plaquette operator for a lattice made of a single plaquette as in Fig. 1 1 v). It 
turns out that 

Z/jZ/ 2 Z/ 3 Z/ 4 : |f,W;) — > |f,-+if/+if/+ifj+i), (B16) 
Z\Z\Z\Z\ A : l^fM-) -> i*-i_i> . (B17) 

It is easy to realize that the eigenvalues of Hp are —1 or +1 /2 and its ground-state on a single plaquette is given by 

\Ya) = -j=Y i \mtiti). (B18) 

V3 ; 

Here, again we discover the emergence of gauge invariant plaquette states. Their number is now 3 and clearly coincide with the 
\G\, the cardinality of the group. The plaquette Hamiltonian is still gapped, but the gap, A3 = 3/2 is now smaller than in the Z2 
case. We can again study the opposite limit of (18) with 9 = n/2. There we see that all spins are aligned in the t\ direction. 
As before, by passing to larger lattices, we can repeat the same reasoning and we find two different phases. One in which the 
ground-state in a product state and the other in which it is robustly entangled. It is also known that by increasing the rank of the 
group a third phase appear. Here, we will not study that phase and refer the interested reader to the literature [67]. 
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